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One of the main tasks for present and future dark energy surveys is to determine whether the 
dark energy is dynamical or not. To iUustrate this from data, it is commonly used to parameterize 
the dark energy equation of state uj as several piecewise constant uiiS using the principal component 
analysis (PCA) method over finite redshift bins. We show that there is only j — 1 free parameters 
uJiS if we choose the redshift as j bins. Without this constrain, one obtains the inconsistent results 
from the data analysis. Furthermore, if uj decreases with non-negligible ratio as z does, then PCA 
fails to reproduce the original behavior of uj. Also, time varying u) can be confused with the incorrect 
f^ . value of constant one when the decreasing (or increasing) ratio of uj is small but not negligible. 
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PACS numbers: 95.36.+X, 95.80.+p, 98.80.Es. 



One of the possible explanations for the current accelerating expansion of the Universe is the existence of the exotic 
energy budget (i.e. dark energy) in addition to the matter component. The different dark energy models are classified 
XV . by the corresponding equation of state (EOS) w, defined by the ratio of its pressure and energy density. Due to 
our ignorance of the nature of dark energy, it might be appropriate to study the dark energy phenomenologically by 
using the accumulating high precision observational data. One needs to parameterize w in order to fit the related 
parameters to data. One general way to parameterize ui is to approximate it using the piecewise constant bins [i , _']. 
C) ] It is claimed that one can reconstruct the time dependence of uj and make further model independent studies by using 
_; ■ a principal component analysis (PCA) method. When we consider the time varying EOS, we adopt the so-called 
f-^ ] Chevallier-Polarski-Linder (CPL) parameterization w = cjq + Waj^ and use its present value (wq) and time variation 
I ■ {uja) [•), 4]. One of the main reason for using this PCA method is known as to determine whether the dark energy 
density is evolving with time or not. However, we show that the time varying uj can be confused with the incorrect 
value of constant one when the decreasing (or increasing) ratio of uj is small but not negligible (i.e. when w^ < 0.3). 
jrt " Thus, even if the PCA gives the similar result as the constant EOS result, we still can not rule out the time varying 
—^t UJ with non-negligible value of Wa model. Also, when the time variation of a; is significant (i.e. when Wq > 0.3), then 
PCA analysis can hardly reproduce the original character of uj. 

We first want to emphasize the correct degree of freedom in PCA. One divides the redshift range of the survey 
f— ^ ■ (z = 0, Zmax) into N bins of not necessarily equal widths Az^ (i = 1, • • • , N), where ^^ Az^ = z^ax- Then, it is well 
f~^ \ known that a set of N values of observations of possibly correlated variables can be orthogonally transformed into 
f — ■ ■ a set of j values of uncorrelated variables so-called principal components [5, 6]. The dark energy is parameterized 
in terms of w(z), which is defined to be constant in each redshift bin, with a value uji in ith bin. For the piecewise 
\f-s . constant uj{z), the energy density of the dark energy for z in bin j evolves as 

o 
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Pde(^) = Pde(^ = 0)^^ _^ii , (1) 



S^ ■ where Zi is the lower redshift bound of the ith bin and uji is the fiducial value of the EOS in that bin. The common 
^ ' mistake in literature is that one regards that there are j free parameters after one fixes z^s. However, we should 
emphasize that aj,;s have only j — 1 degree of freedom because pde{z) should be equal to pde{z = 0) when z = 0. 
From the above equation (1), this is given by 







l+(^(l+...)ln[i±f|i|)/(ln[l + z,]). (2) 



Thus, ujj is determined by other parameters uji, Zi, and Zj. Of course, one can constrain one of w^s instead of ujj 
without changing the final result. 
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TABLE I: Parameter values of true models are shown in the first two columns. Zi is the uncorrelated bin and uJi and uj* are 
obtained from the minimum y^ fitting by using the simulated H{z) data compared to the true model with and without the 



constraint on the one of parameter values, respectively, 
values with and without the constraint. 



(Ji is the 1-(J error of uji. x a-nd x* correspond to the minimum x 



In what follows, we use the 40 simulated equally binned Hubble parameter H{z) data for < ^ < 2 by using the 
true models given by CPL parametrization. We assume that the measurement error on iJ as 5 % and there are no 
errors in both the present energy density contrast of the matter fij^ and the present Hubble parameter value Hq . We 
also consider only the flat universe. We perform a simple xlf test to determine the best fit values of uJi and \-a error 
of each Wi is obtained from the covariance matrix. We check the reliability of PCA method to constrain the DE and 
using H{z) is good enough for this purpose. When we investigate the luminosity distance, the results are even worse 
because of the multi-integral of w in it. We show the results for the two different models in table I. In the first case, 
true model is characterized by (wo,Wa) = (—0.9,0.5). We perform y^ test with H{z) data created from the this true 
model with the assumption of 5 % error. The true model of the second case is given by (wcWq) = (—1.1,0.3). We 
explain the details of the data with Fig. 2 later. There is one remark to be emphasized. The 1-cr error increases as Zi 
does. This is due to the fact that x^ fitting is performed from the lowest Zi. Thus, the errors in lower z^s is propagated 
to the higher ZiS. Thus, no matter how much data point we add in the higher z^, this intrinsic error propagation will 
not be removed. We choose the uncorrelated bins as z = 0.1, 0.4, 0.8, 1.25, and 1.6. As we show before, there are 
only 4 degree of freedom in this case. Thus, the value of W5 at z = 1.6 is derived from the other uJi and Zi values by 
using Eq. (2). We compare the values of w^ with the ones of uj* which is obtained from x^ test without the above 
constraint. As we can see there are discrepancies in the uj^ values between them. 

Even though, the above constrain Eq. (2) seems to be trivial but the piecewise parametrization of a; suffers the 
inconsistency without this. For example, one obtains the best fit values of tUi from the x^ fitting after z^, Zj, il^, and 
-ffo are specified. Without the above constrain Eq. (2), one is not able to recover the original 51^ value which is used 
in the maximum likelihood analysis. This is true whether one uses the marginalization of nuisance parameter 17^0 
in the analysis or not. We show this in Fig. 1 by using w^ and zi values in table I. When one use the raw obtained 
value of w* (dotted lines), ri°j is different from the one used in the analysis as shown in the Fig. 1. The evolutions 
of Q.m{z)s of the true models are depicted as the dashed lines. The evolutions of D,„i{z)s obtained from the correct 
consideration of degree of freedom are shown as the solid lines. The left and right panels correspond to (wq, w„) = 
(—0.9,0.5) and (—1.1,0.3), respectively. 



We find several drawbacks in PCA. Firstly, when the time variation of EOS is not negligible, PCA can not produce 
the proper behavior of EOS at the entire region of the investigated rcdshifts. Secondly, PCA produces the EOS which 
is confused with the constant EOS with the improper present value of it when it changes slowly (i.e. Wq < 0.3). In 
Fig. 2. we demonstrate both cases. In the left panel of Fig. 2, we show llh (solid) and w* (dotted) when the true model 
(dashed) is (a;o,Wa) ~ (—0.9,0.5). The obtained values of Wi show the oscillation behavior around —0.7 even though 
the true model decreases monotonically. Thus, PCA method produces the totally different behavior of w when Wa is 
non-negligible. Also in the right panel of Fig. 2, the obtained w, is almost same as that of the cosmological constant 
(A) for the entire region of z even though the true model is (wq, Wa) — (~1-1, 0.3). We can compare this result with one 
in Ref. [7]. Even though the result in the mentioned reference seems to be consistent with the cosmological constant, 
there still can be the viable time varying DE models which can mimic A. This impedes any proper interpretation of 
the result obtained from PCA method. Thus, PCA also may mislead to the true property of dark energy in the slow 
changing a;. We check that PCA method can give the reliable result only when w is almost constant. 





FIG. 1: The evolutions of Qrn{z) obtained from a;* (dotted), ljo + ^ajir^ (dashed), and Ui (sohd). a) When {uojUJa 
(-0.9,0.5). b) For (^0,^^^) = (-1.1,0.3). 
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FIG. 2: Comparison between a;, oji, and oj*. a) True model is w = —0.9 + 0.5^4^ (dashed) and the obtained values of lOi from 
PCA with the constraint (solid) and without the constraint on the last component (dotted), b) True model is ij; = — 1.1 + 0. Syr— 
(dashed) and the obtained values of a;^ from PCA with the constraint (solid) and without the constraint on the last component 
(dotted). 

Even though, PCA is the most model independent method for probing DE, we show that the true degree of freedom 
of the parameters should be one less than the binned number. We also demonstrate that PCA method may mislead to 
the property of dark energy when the time variation of it is not negligible. Even for the slowly varying w, PCA result 
may produce the incorrect information on the true DE. We check that the above conclusion is same for the increasing 
a;. PCA is adequate only when a; is a constant. Thus, we may need to check both CPL like model dependent w 
parametrization and PCA method to investigate the DE properly. 
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